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ABSTRACT 

We consider envelopes of DB white dwarfs which are not composed of pure "^He, but rather a mixture of ■'He 
and "^He. Given this assumption, the same diffusive processes which produce a relatively pure H layer overlying 
a He layer in the DA's should work to produce a ^He layer overlying a "^He layer in the DB's. We examine the 
relevant timescales for diffusion in these objects, and compare them to the relevant evolutionary timescales in the 
context of the DBV white dwarfs. We then explore the consequences which ^He separation has on the pulsational 
spectra of DBV models. Since GD 358 is the best-studied member of this class of variables, we examine fits to 
its observed pulsation spectrum. We find that the inclusion of a ^He layer results in a modest improvement in a 
direct fit to the periods, while a fit to the period spacings is significantly improved. 

Subject headings: dense matter — stars: oscillations, evolution — white dwarfs 



1. ASTROPHYSICAL CONTEXT 

The two greatest successes of white dwarf asteroseismology 
have come from the analysis and interpretation of Whole Earth 
Telescope (WET) data on the objects GD 358 and PG 1159 
(Winget etal. 1991, 1994). The temporal spectra of both objects 
show well-defined multiplet structures as well as many consec- 
utive radial orders of the same t. This allows us immediately to 
estimate the rotational frequency, VL, and the mean period spac- 
ing, {AP)(, with the latter quantity putting strong constraints 
on the mass of the best-fitting model. Furthermore, the vari- 
ations from equidistant period spacing, i.e., "mode trapping", 
give us information (at least in the case of GD 358) about the 
radial structure of the star, such as the thickness of its surface 
He layer. 

In a detailed analysis of the frequency distribution in GD 358, 
Bradley & Winget (1994) found evidence of a chemical tran- 
sition zone located at a mass depth of 1.5 x 10~^M,t, which 
they interpreted as the C/He boundary. A He layer thickness 
of ~IO~^M* is thinner than expected on naive evolutionary 
grounds, as well as thinner than suggested for the DA white 
dwarfs by the asteroseismological work of Clemens (1993). 
While the fit of Bradley & Winget explained the trapping fea- 
tures in the neighborhood of 700 sec reasonably well, it did not 
do as well for periods near 500 sec. They found evidence that 
these periods near 500 sec might be better fit if an additional 
chemical transition zone at a depth of ~10~^M,t is assumed. 

Independent of the above considerations, if the theory of dif- 
fusion of chemical species in stars is correct, then at least a par- 
tial separation of different isotopes of the same element must 
occur in white dwarfs, as was pointed out to one of us by Clay- 
ton (1988, private communication). At the time, there were no 
known observational consequences; this has since changed with 
the advent of white dwarf asteroseismology. 

While the Galactic number ratio of ^He to ''He has been mea- 
sured to be of order lO"'' in many astrophysical environments 
(Galli et al. 1995; Prantzos 1996), this does not necessarily hold 



for the stellar cores which are presumably the white dwarf pro- 
genitors. In fact, standard evolutionary theory would suggest 
quite low levels of ^He in these objects (Galli et al. 1995). Our 
view is that given the theoretical uncertainties in the late stages 
of stellar evolution, we should not rule out any possibilities, 
and indeed should seek to make as many independent measure- 
ments as possible. 

These considerations led to two independent motivations for 
considering ''He diffusion (Winget 1998). First, it is a process 
which may be a generic feature of white dwarf cooling, and, as 
such, should be considered as part of our asteroseismological 
analyses. Second, it might allow us to place the mode trapping 
data for GD 358 in a different context: this star could have a 
C/'*He transition zone at ^ lO'^M,^ and a thinner ''He/^He tran- 
sition zone at 10~*M*. If this is the case, then the DBV's 
and the DAV's would again have the same order of magnitude 
He layer thicknesses. To determine the plausibility of these hy- 
potheses, we first examine the relevant diffusion timescales, and 
then the effect which such a layering structure has on our fits of 
GD 358. 



2. DIFFUSION 
2.1. Timescales 

The process of diffusion is of major importance, since we 
are aware of no other effect which could lead to the spatial 
separation of ^He and ''He. We therefore examine the relevant 
timescales for such diffusive processes, to see if it is plausible 
for a significant fraction of the ^He to have separated from ''He 
in the elapsed evolutionary times for these objects. 

Fortunately for us, Fontaine & Michaud (1979) have already 
examined the related problem of C diffusion in a background 
of normal He, i.e., ''He. In their analysis, they treated C as a 
trace element. This is an excellent approximation for our case 
as well, since we naively expect that the ^He is only about one 
part in 10'' of the "He. 
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Fig. 1 . — The dmescale, 9, for diffusion of ^He in a background of ''He as a 
function of T^ff at different mass depths between -12 and -2. The curve labeled 
"Age" gives the age of the model as a function of Teff; belovif the line the profile 
should have reached diffusive equilibrium and above it little diffusion should 
have occurred. 



3. EQUATION OF STATE 

To include the effects of a ^He layer we have made a rela- 
tively simple modification to the envelope routines in our evo- 
lutionary code. Since we are interested in modehng a DB, we 
are free to use the array space normally reserved for the H pro- 
file and use it for the ^He profile. To this end, we have replaced 
the H equation of state (EOS) and opacities with those appro- 
priate for ^He. In doing this, we have taken the ^He EOS and 
opacities to be equal to those of "^He at 3/4 the density (to cor- 
rect for the isotopic mass ratio). This approximation should be 
more than sufficient for our purposes. A comparison of the re- 
gion of period formation for the cases of a "^He/H and a "^He/^He 
envelope is given in Montgomery & Winget (1999). 

4. PREVIOUS BEST FITS 

We now wish to examine the effect which a ^He layer could 
have on the pulsation frequencies. First, in the top panel of Fig- 
ure 2, we show one of the best-fit models for the star GD 358 
from Bradley & Winget (1994). The filled circles connected by 
solid lines show the observed mode trapping structure, and the 
open circles connected by dotted lines show the results from the 
best-fit model. This model has M^/M© = 0.61, T^ff = 24,044 
K, and M4He/M^ = 1.5 • 10"*'. The lower 3 panels show the 
effect a thin layer of ^He has on the mode trapping structure. 
If the layer is as thick as M^ue/Mj, = 1.5 • 10"'", corresponding 
to A'(^He)/A'C*He) = 10"^, then we see that the mode trap- 
ping structure in the vicinity of 700 sec is significantly altered. 



Applying their equations for this case for a 0.612 Mq model, 
we obtain the result shown in Figure 1. We see that if we as- 
sume a Teff of 25,000 K for a typical DBV, then diffusive equi- 
librium between ^He and "^He should prevail down to approxi- 
mately the lO^^M* mass point. Thus, if the "^He layer is as thick 
or thicker than this, and the ^He abundance is ~ 10~^ that of 
the ''He abundance, we would expect a ^He layer of thickness 

~ lo-l 

One clear prediction of this is that the spectra of DB's in 
this temperature range should show more or less pure ^He, at 
least for ''He layers which are thicker than ^ lO"*". This is cur- 
rently in the process of being tested using high-resolution spec- 
troscopic observations (Koester, private communication). 

We mention as an aside that the case of ^He may be unique 
in that it is a trace isotope which is lighter than the dominant 
species, so that it can diffuse upwards and produce a thin, but 
not asteroseismologically negligible layer. For trace isotopes 
which sink, there should be no such signature. 



2.2. Profiles 

In this preUminary exploration, we will only consider 
'*He/^He zones with equilibrium profiles, although the profiles 
could in fact be less sharp than this. The treatment of the chem- 
ical profiles in the transition zones in our present and previ- 
ous work is based on the work of Arcoragi & Fontaine (1980). 
Essentially, we use equation (A6) of Arcoragi & Fontaine 
(1980), which assumes that an equilibrium distribution has been 
reached and that one of the elements may be treated as a trace 
element. This is certainly a valid assumption for the initial 
stages of ^He diffusion, given the expected abundance ratio of 
^He to ''He. 
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Fig. 2. — A comparison of the observed trapping (filled circles, solid lines) 
with carbon core models (open circles, dotted lines). The top panel is the best- 
fit model of Bradley & Winget (1994), which of course has M^ng = 0.0. The 
lower panels show the mode trapping structure for nonzero values of Mjue, as 
indicated in each panel. 
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Table 1 

Results of Fits using a Genetic Algorithm 



Model 


Tiff (K) 


M^/Mq 


l0gM4He/M* 


logM3He/M* 


aiP) 


a(AP) 


C/thick4He 


23,100 


0.655 


-2.74 




2.30 


2.65 


C/thick4He/3He 


23,200 


0.605 


-2.74 


-6.02 


1.99 


2.29 


C/thin4He 


23,800 


0.600 


-5.76 




2.71 


2.53 


C/thin4He/3He 


24,400 


0.580 


-5.71 


-10.09 


1.95 


2.66 


0/a4He'' 


22,600 


0.650 


-2.74 




1.50 


1.92 


0/a4He/3He'' 


22,300 


0.630 


-2.79 


-7.49 


1.29 


1.32 


Inferred Errors 


500 


0.015 


0.06 


0.12 







"These fits were performed using / = 1 / [u{P) + ct( AP)] as the fitness criterion; for the other models f -1/ <t{P) was used. 



Unfortunately, this is the region in which the fit was already 
quite good, whereas in the region around 500 sec, where the 
fit was poor, there is very little change. However, this analysis 
clearly demonstrates that the inclusion of a ^He layer can have 
a measurable effect on the calculated mode trapping structure. 

5. NEW FITS USING A GENETIC ALGORITHM 

In an effort to test asteroseismologically the viability of the 
hypothesis of a ^He layer, we performed an extensive set of cal- 
culations. Using a genetic algorithm (Metcalfe 1999; Metcalfe, 
Nather, & Winget 2000), the parameters Mi,, Teff, M4He, and 
M3He were varied to produce pulsational models whose l=\, 
m = Q periods were compared with those inferred from the ob- 
servations. The region of parameter space explored is given by: 
0.45 M0 < < 0.95 Mq 

20,000 K < Teff < 30,000 K 
~10-^M^ < M4He < 10-2 M,^ 

10-5 M4He < M3He < 10-^ M4He 

The best-fit models from these runs are listed in Table 1 . The 
column labeled aiP) gives the residuals (standard deviation) of 
the calculated and observed periods, and the colunm labeled 
a{AP) gives the residuals of the calculated and observed period 
spacings. For the upper four models we have taken the fitness 
criterion to be 1 /(t{P), whereas for the lower two models we 
took it to be 1 /[cr(P) + fT(A/')]; thus, for these last two entries, 
we are fitting not just the periods but the spacings between con- 
secutive periods. We have chosen this criterion for the models 
with a more realistic core composition based on our suspicion 
that the period spacings will be a more sensitive diagnostic of 
the ^He layer than the periods themselves will be. 

As a reference to previous fits, the fit of Bradley & Winget 
(1994) for a C core model has a(P) ^ 2.3 sec. Thus, the present 
best fits, both with and without a ^He layer, having period resid- 
uals of ~ 1.5 and 1.3 sec, respectively, represent a significant 
improvement over previous fits. 

6. DISCUSSION 

We see from Table 1 that the fits with O/C cores are signifi- 
cantly better than those with pure carbon cores, independent of 
whether a ^He layer is present or not, a result previously found 
by Metcalfe, Nather, & Winget (2000). In Figure 3, we display 
the mode trapping diagrams for these two fits. Both fits repro- 
duce the periods quite well, with the ^He model reproducing the 



period spacings much better (~ 1 .3 sec compared to ~ 1 .9 sec). 
This indicates that the period spacing may be a better diagnostic 
for the fine structure produced by the "^He/^He transition zone 
than just the periods themselves. 

We now seek to understand the relative importance of the 
O/C and "'He/^He transition zones. In the asymptotic limit of 
high radial overtones and large periods, the frequency of a given 
mode is given by a simple radial integral of the Brunt-Vaisala 
frequency. Using a "period formation" diagram, we can show 
the relative weight which a given region has in determining a 
mode's period. 

In Figure 4, we show such a diagram for the case of our 
best-fit ^He model. The three peaks which are labeled corre- 
spond to the O/C, C/'*He, and ''He/^He transition zones. Using 
the equilibrium diffusion coefficients for the "^He/^He transition 
zone, we see that it is of relatively minor importance in deter- 
mining the mode frequencies. In contrast, the O/C and C/*He 
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Fig. 3. — The upper panel shows the best-fit model without ^He, and the 
lower panel the best-fit model with 'He. As can be seen, these fits have much 
smaller residuals than the previous fit given in the upper panel of Figure 2. 
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transition zones are both quite pronounced and should signifi- 
cantly affect the periods of the modes calculated in the models; 
this is borne out by the major improvement in the standard 
deviation of the periods, a{P), with the inclusion of an O/C 
chemical profile. The period spacings, on the other hand, are a 
differential quantity, and therefore sensitive to even small devi- 
ations in the background structure. It is therefore not surprising 
that the inclusion of a ^He layer results in a major reduction in 
the resduals of the period spacings. 

We now attempt to quantify the statistical significance of the 
improvement of the fit when the parameter corresponding to the 
^He layer thickness is added. Following Koen & Laney (2000), 
we apply the Bayes Information Criterion (BIC) for = 1 1 data 
points. We find that the addition of a parameter should be ac- 
companied by a decrease in the residuals of at least ~10% in 
order to be considered statistically significant. We note that this 



is a necessary but certainly not sufficient condition for the va- 
lidity of adding a parameter. 

Examining the best-fit O/C core models, we see that the stan- 
dard deviation of the fit to the period spacings improved from 
1 .92 sec to 1.32 sec with the addition of ^He. This is a decrease 
of 30% and should therefore be considered statistically signifi- 
cant. 

Finally, we again note that the models with O/C cores yield 
markedly lower residuals, regardless of which criterion is used 
(periods or period spacings). This result has recently been ob- 
tained by Metcalfe, Nather, & Winget (2000), who report the 
results of such fits in an extensive parameter space of white 
dwarf models. 

7. CONCLUSIONS 

From this preliminary analysis, we have shown that diffusion 
theory applied to white dwarfs predicts that any ^He initially 
present in the DBV's down to a depth of ~10~^Mi, should have 
diffused upward to produce a surface layer of ^He. In addition 
to being detectable spectroscopically (Koester, private commu- 
nication), such a ^He layer can significantly affect the astero- 
seismological fits, and therefore needs to be included in such 
analyses. 

We find that although the inclusion of a ^ He layer results in 
only a marginal improvement to the fits to the periods, the fit 
to the period spacings is significantly improved. This is be- 
cause the period spacings are more sensitive to the fine struc- 
ture which a ^He/*He transition zone produces than the periods 
themselves are. 

Finally, we find that an O/C core (essentially, the transition in 
the core from an O/C mixture to pure C) fits the observed pul- 
sational spectrum of GD 358 much better than a pure C core, in 
agreement with Metcalfe, Nather, & Winget (2000), who first 
obtained this result. This gives us the hope of someday being 
able to constrain the prior nuclear burning history of GD 358 
and other pulsating white dwarfs. 
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